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O . Abstract 

Q^ ■ It is shown that the fine structure sphttings of the 2 Pj and 3 Pj 

excited states in charmonium are as large as those of the l^Pj state 
CL|| if the same as(^) ~ 0.36 is used. The predicted mass M{2^Pq) = 

p ■ 3.84 GeV appears to be 120 MeV lower that the center of gravity of 

D , the 2^Pj multiplet and hes below the DD* threshold. Our value of 

"^ ' M{2^Pq) is approximately 80 MeV lower than that from the paper 

by Godfrey and Isgur |1| while the differences in the other masses 
^ . are < 20 MeV. Relativistic kinematics plays an important role in our 

H I analysis. 

1 Introduction 

At present only the l^Di and 2 ^Di states lying above the DD threshold 
have been identified with the experimentally observed cc mesons, iIj{3770) 
and iIj{4160). Still a large number of other excited P- and D-wave states 
above the fiavor threshold were predicted. Their masses and fine structure 
splittings were calculated by Godfrey and Isgur (GI) already in 1985, in the 
framework of a relativistic approach IQ. The properties of the P and D 
levels in charmonium and bottomonium were also intensively studied in the 
nonrelativistic approximation [^, ^. There is a point of view that one or 
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more charmonium 2Pj states can be narrow enough to have a substantial 
branching ratio to the 7 + 'ip{2S) channel ||^, ^] and could play a role in the 
hadronic production of ip{2S) mesons. In particular, the 2P states can be 
related to the enhancement in the J/t/'tt+tt" system near M = 3.84 GeV 
observed in P] (but not confirmed by another group 0). Therefore the 
precise knowledge of their masses is especially important. 

An accurate description of the charmonium spectrum and the fine struc- 
ture splittings of the IP level was presented in our previous papers p, |[, 
where, as well as in [|l|, relativistic kinematics was taken into account with 
the help of the spinless Salpeter equation. As shown in |0 the relativistic 
corrections to the matrix elements, like < r~^ > determining spin structure, 
are large enough, of the order of ~ 40%, therefore the nonrelativistic ap- 
proach can not be considered to be appropriate when the spin structure is 
investigated. 

We shall show here that in the relativistic approach the spin-orbit and 
tensor splittings are rather large for the P states, for the ground state as well 
as the excited levels. This result depends weakly on the choice of the strong 
coupling constant ^^(/i). Here as(/i) ~ 0.36 (/i = 0.92 GeV) will be used for 
all states, but the splittings do practically not change if as{fi) = 0.30 (/i = 
m = 1.48 GeV) is taken. 

The fine structure splittings predicted here appear to be larger than those 
in GI's paper |l|], especially for the 2P and 3P states. The reasons for that 
will be discussed in Section 3. As a result the 2 ^Pq mass, M(2 ^Pq) = 3.84 
GeV in our case appears to be approximately 80 MeV smaller than in |jl| 
and this level lies below the DD* threshold. The 2 ^Pi,2 levels as well as the 
n ^Dj states have mass values close to the GI predictions. For the first time 
we also predict large fine structure splittings for the 3P states. 

2 Spin— averaged spectrum 

The relativistic effects in charmonium are expected not to be small, especially 
for the wave functions and matrix elements which we are mostly interested 
in here. Therefore to find the spin-averaged spectrum the spinless Salpeter 
equation will be solved as already done in several papers |l|, [1^, |Tl[ 

(^2^p^ + m^ + Vo{r)^ ^ni{r) = M^i^^i{r). (1) 



The static interaction Vo(r) will be taken in the form of the Cornell potential, 

Voir) = -^- + ar + Co (2) 

3 r 

and the values of a = ay(/i), the string tension a and the pole mass of the c 
quark will be taken as in our previous paper p while fitting the fine structure 
of the IP charmonium state, 

m = 1.48 GeV, a = 0.18 GeV^ a = 0.42. (3) 

The constant Cq in Eq. (2) is determined from a fit to the spin-averaged mass 



of the 15 states, M{IS) = 3067.6 MeV p[, and from there Co = -140.2 



MeV. In our approach the strong coupling constant a is not an independent 
parameter. It can be connected to as{fi) (in the MS renormalization scheme) 



through the relation [1^ 



a{q'^) = aviq'^) = asifi) (l + - — as(/u)j (% = 3). (4) 

For our choice of a and ^^(/i) (see below) this relation will be valid with an 
accuracy better than 5%. 

The chosen parameters (3) can be compared to the GI parameters in 
||l]| where m = 1.628 GeV while a = 0.18 GeV^ coincides with a in (3). 
In |T]] a running coupling constant was used for a with the critical value 
acr = cy{r = 0) = 0.60 which is larger than the constant a = 0.42 in 
our case. Also in [^ Cq = —253 MeV whereas Cq = —140 MeV in our 
calculations. Nevertheless, the calculated spin-averaged masses for both sets 
of parameters are close to each other. For P-wave states the differences are 
less than 10 MeV and for D-wave states less than 20 MeV (see Table |l|). 

In many papers excited states of charmonium were studied in the non- 
relativistic approximation which works quite well for the spectrum. But it 
was shown in P, ^ that the relativistic corrections to matrix elements like 
< r~^ > and < r~^ In mr >, which determine fine structure splittings, are 
rather large, about 30-i- 40%. That is why in this paper only relativistic cal- 
culations of the fine structure splittings of the excited states of charmonium 
will be considered. 



3 Fine structure parameters of the P levels 

Although the spin-averaged masses in our calculations are very close to those 
in GI's paper [|l|], we expect that the spin-orbit and tensor splittings of the 
P-wave states will be larger in our case. There are two reasons for that. 
First, we take into account the corrections to second order in a^. Secondly, 
our calculations of different matrix elements have shown that for the excited 
P states the matrix element < r~^ > which determines the splittings, is not 
decreasing. For the set of parameters (3) it was found that in the relativistic 
case < r-'^ >ip= 0.142, < r"^ >2P= 0.157, < r~^ >3P= 0.167, i.e. < 
r~^ >nP are even increasing for 2P and 3P states. This result is specific to 
the Salpeter equation, whereas in the nonrelativistic case the matrix elements 
< r~^ > for the excited states are decreasing, e.g. < r~^ >ip= 0.101, < 
r"^ >2P= 0.093, < r~^ >3P= 0.089. The accuracy of our calculations was 
checked to be (1^2)10-1 

The fine structure parameters are defined as matrix elements of the spin- 
orbit and tensor interactions, 

a =< VLsir) >, c =< Mr) >, (5) 

where the scalar functions VLs{r) and Vr(r) are introduced as 

VLs{r) = VLs{r)L-S, 

~ « « r 

Vrir) = VT{r)Si2, Su = ^{s- n){s2 ■ n) - si ■ S2, n = -. (6) 

r 

Here the spin-orbit parameter a is defined in the same manner as in other 

papers, whereas the definition of the tensor parameter c differs from |[l|] where 

the tensor parameter T = ^c and from where the parameter b is used, 

related to c by 6 = 4c. 

In our calculations we suggest that the P— wave hyperfine splitting is 

small as it occurs for the he (IP) meson for which the hyperfine shift relative 

to the center of gravity of the ^Pj multiplet, M(l'^Pj), is less than 1 MeV. 

When the hyperfine splitting is neglected the mass of the states with spin 

S = coincides with the center of gravity of the ^Lj multiplet denoted by 

Ml. Their values, taken from Table |I], are 

M(l^Pi) = 3528 MeV; M(2 ^Pi) = 3962 MeV 

M{1^D2) = 3822 MeV; M(2 ^Ds) = 4194 MeV, (7) 



which are about 20 MeV lower than in |T| for the D and some S states. 

For the states with spin S = 1 and orbital angular momentum L 7^ the 
mass of the state can be represented as 

M{^Lj) = ML + a<L-S> +c<Su> (8) 

with the operator 5*12 defined as in (6). For P-wave states it gives 

M{^P2)=Mi + a-—c, M{^Pi)=Mi-a+-c, M{^Po) = Mi-2a-c. 

(9) 
In Ref. |1| only the terms of first order in ag were taken into account and 
for a and c the following values were obtained, 

a(lP) = 28 MeV, c(lP) = 26 MeV, 

a(2P) ^ 17 MeV, c(2P) ^ 8 MeV, (10) 

which are about 20% and 30% resp. smaller for the IP state compared to 
the current experimental values [§: 

aexp(l^) =34.56 ±0.19 MeV Cexp(lP) = 39.12 ± 0.62 MeV. (11) 

In our approach the terms of second order in a^ will be taken into account 
and the total value of a and c can be represented as 

fltot = Cip^ + flp^ + «NP, Ctot = C'p + c'-p + Cnp, (12) 

where the nonperturbative contribution to the spin-orbit splitting coming 
from the linear confining potential is 

a,p = -£,<r-'>- (13) 



In the tensor splitting, Eq. (\l2\j, the small nonperturbative term Cnp will be 
neglected, (see the discussion in ||^). Godfrey and Isgur took into account 
non-perturbative spin effects. In order to do so, they needed to smear their 
potentials at short range, which makes it necessary to intoduce additional 
unknown parameters. Here we consider spin-effects as a perturbation using 
explicit analytical expressions for the spin-orbit and tensor potentials in 



coordinate space in the MS renormalization scheme from []T^. The terms to 
first order in a^ are Op and Cp given by 

a« = ^<r-3>, c« = ^^^ (14) 

m^ 3 m^ 

and the second order perturbative corrections are 

f2) ^ a^[ii) r^^^ 1j^ ii < ^-3 > _^2.5 < r-3 In(mr) > +1.582 < r"^ >1 



^(^)-^{4.51nii/^-3 



c^^ = _£^{4.5ln£ < r""* > +1.5 < r"Mn(mr) > +3.449 < r'-^ >} 

(15) 

The second order expressions are given here for a number of fiavors nj = ?>. 
(We checked that for n/ = 4 the values of a and c do not practically change- 
the differences are less than 0.5 MeV-therefore only the case nj = 3 will be 
presented here.) 

With the solutions of the Salpeter equation (1) all matrix elements deter- 
mining Eqs ( p!3| - p!5| ) can be calculated and the only uncertainty comes from 
the choice of the strong coupling constant ^^(/x) and the value of the renor- 
malization scale /i. In |^ it was found that for the IP state in charmonium 

a,(/i) = 0.365 (/i = 0.92 GeV) (16) 

gives an accurate description of the spin splittings. For the excited 2P and 
3P states where experimental data are absent, we shall use the same value 



([T6| ) for as{fi)- The main argument in favour of this choice can be taken from 
the fine structure analysis in bottomonium where the values of ^^(/i) for the 
IP and 2P states differ by about 10% only M. 



With asdj-) = 0.365 the authors of |]T4[ obtained for the IP state 

a5,^)(lP) = 47.6 MeV, ag^(lP) = 3.6 MeV, aNp(lP) = -16.6 MeV, 

c^p\lP) = 31.7 MeV, cg^(lP) = 7.4 MeV, (17) 

so that atot(lP) and Ctot(lP) just coincide with their experimental values 



m 



For the excited 2P state the spin-orbit and tensor parameters are as 
large as those of the IP state, because the matrix element < r~^ >2p is 
even ^ 10% larger than < r~^ >ip- Here we face the difference between the 
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relativistic and nonrelativistic approaches. For the latter, matrix elements 
like < r~^ >np are decreasing with growing n = n^ + 1. Our calculations 
give for the 2P state that 

af{2P) = 52.5 MeV, ag^(2P) = -0.4 MeV, a^p{2P) = -13.4 MeV, 

c'^p\2P) = 35.0 MeV, c'-p\2P) = 6.5 MeV (18) 

so that 

atot(2P) = 38.7 MeV, Ctot(2P) = 41.5 MeV (19) 

are even slightly larger than the corresponding values for the IP state. 

Comparing values obtained for a(2P) and c(2P) with those in (11) from 
GI's paper one can see that a and c in our calculations are larger by a factor 
of 2 or 5 resp. than the values from GI's paper (see also Table 0). This 
discrepancy is partly connected to our taking into account of the second 
order radiative corrections which are, however, not large. But even with 
only first order perturbative terms included, our values of a(2P) and c(2P) 
are much larger than GI's values (11). 

Taking a and c from Eq. (|19|) and Mi(2P) = 3962 MeV (from Table D 
one obtains the following masses for the 2 ^Pj states, 

M(2 ^Po) = 3843 MeV, M(2 ^Pi) = 3944 MeV, M(2 ^Pa) = 3997 MeV. 

(20) 
Our predicted mass for the 2'^Po state, M(2^Po) = 3.84 GeV appeared to 
be ~ 80 MeV lower than the one in GI's paper, whereas for the other two 
states, 2^Pi and 2'^P2, the predicted masses only slightly differ from the GI 
values (see Table 0) due to a cancelation of terms with different signs. 

It is important that in our calculation the 2 ^Pq level lies below the DD* 
threshold {Mth{DD*) ^ 3.87 GeV) but higher than Mth{DD) = 3.73 GeV. 
This fact can affect the decay rates of the 2 ^Pq state. 

For the 3P state, using again ^^(/i) = 0.365 and fi = 0.92 GeV as for the 
IP state, one obtains using Eqs. (|T3| - |T5D 



a5,^)(3P) = 56.7 MeV, ag^(3P) = -3.6 MeV, aNp(3P) = -11.8 MeV, 

4^) (3P) = 37.8 MeV, cg^(3P) = 6.2 MeV, (cnp = 0), (21) 

so that 

a,^,(3P) = 42.3 MeV, Ctot(3P) = 44.0 MeV. (22) 



With these values for a and c and the spin-averaged mass Mi(3P) = 4320 
MeV it follows that 

M(3 3Pi) = 4300 MeV, M(3 ^Pa) = 4358 MeV, M(3 ^Pq) = 4192 MeV. 

(23) 
The level 3 ^Pq lies 128 MeV lower than the center of gravity of the 3 ^Pj 
multiplet. It is of interest to notice that the difference M{n '^Pa) — M(n^Po) = 
A{nP) = 3a + 0.9c is large in all cases and slightly increasing for the excited 
states, 

A(1P) = 138.9 MeV, A(2P) = 143MeV, A(3P) = 166MeV. (24) 

We have checked the sensitivity of the predicted values for a and c to the 
choice of the renormalization scale fi and as (fJ^) ■ To this end the commonly 
used value of /i = m and as{^ = m = 1.48 GeV) = 0.29 was considered. 
Then for the 2P state atot(2P) = 36.5 MeV and Ctot(2P) = 37.5 MeV and for 
the 3P state atot(3P) = 40.5 MeV and Ctot(3P) = 39.8 MeV were obtained, 
which are very close to the values (20), (21), (23), and (24) with /xq = 0.92 
GeV and as{fio) = 0.365, found in §] from the fit to the IP fine structure 
splittings. 



4 Fine structure splittings of the D levels 

For the D state the expressions of the splitted masses, M{n^Dj), through 
the parameters a and c can be found in |^, 

M( ^D^) = M2-3a--c, Mi^D2) = M2-a+-c, M^D^) = M2+2a--c. 

_ (25) 

For the spin-averaged masses M{nD) our calculation with the parameters 

(3) give 

M2(1P') = 3822 MeV, M2{2D) = 4194 MeV. (26) 

The fine structure parameters a and c for the D-wave levels are given in 
Table |] together with their values from GI's paper. As seen from Table |^, 
a and c in both cases practically coincide. Still our predicted masses for the 
I'^Dj and 2 ^Dj states appear to be about 20 MeV lower than in |1[] because 
of the smaller value of the spin-averaged masses. 



5 Conclusion 

In our analysis it was found that 

(i) in the relativistic case the fine structure sphttings of the charmonium 
P-states with spin S = 1 are much larger than in the nonrelativistic case; 

(ii) for the excited 2P states the parameters of the fine structure are even 
slightly larger than for the l^Pj ground state; 

(iii) the mass of the n^Po states {n = 1,2,3) appears to be about 130 
MeV smaller than the center of gravity of the n^Pj multiplet. This fact 
can be important for the explanation of the decays of the excited states of 
charmonium. 

The value we predict for M{2''^Po) = 3.84 GeV is about 80 MeV lower 
than in GI's calculations []I|. This state lies below the DD* threshold and 
is only about 100 MeV higher than the DD threshold. The point of view 
exists that this state could be very broad because it lies above the DD 
threshold and therefore should have a large hadronic width. On the other 
hand this state lies relatively close to the DD threshold, so the width could 
be suppressed by phase space limitations. Therefore this state could play a 
role in the production of ip{2S) charmonium mesons as it was discussed in 
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Table 1: The spin-average masses M{nL) (in MeV) in charmonium for two 
sets of parameters. 



M{nL) 


Godfrey and Isgur 

m = 1.628GeV 

a = O.lSGeV^ 

c^crit = 0.6 (running a{r)) 

Co = -253 MeV 


present paper 
m = 1.48 GeV 
a = 0.18GeV2 
& = 0.42, 
Co = -140 MeV 


Experiment 


^{IS) 


3067.5 


3067.6 


3067.6±0.6 


ij{2S) 


3665 


3659") 


3663±1.3 


ij{3S) 


4090 


4077'') 


4040±10 


^{4S) 


4450 


4425 


4415±6 


^{5S) 




4732 




XcilP) 


3520 


3528 


3525.5±0.4 


Xc{2P) 


3960 


3962 




Xc{3P) 




4320 




M{ID) 


3840 


3822 


3768.9±2.5 


M{2D) 


4210 


4194 


4159 ±20 


M{3D) 


4520 


4519.5 





"•^ Mixing of 25* and ID-wave states is not taken into account 
'') Mixing of 35* and 2D-wave states is not taken into account 
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Table 2: Spin-orbit and tensor splittings a and c (in MeV) for P and D levels. 





Godfrey and Isgur 
paper") 


present paper 
as{^l) = 0.365 


Experiment 


a(lP) 


28 


34.56 


34.56 ±0.19 


c(2P) 


26 


39.12 


39.12±0.62 


a(2P) 


17 


38.7 




c(2P) 


8 


41.5 




a(3P) 




42.3 




c(3P) 




44.0 




a(lP)) 


^5 


3.64 




C(1D) 


^ 10 


10.94 




a(2P)) 


^5 


5.43 




c{2D) 


^ 10 


11.37 





"•'' the values of a, c for 2P and D states are extracted from the masses M{^Dj) 
and M{2^Pj) given in Ref. M. 



Table 3: Masses of the n^Pj, and rr'Dj states (in MeV) in charmonium. 





Godfrey and Isgur 


present paper 






paper 


asifi) = 0.365 


Experiment 


2^Po 


3920 


3843 




23Pi 


3950 


3944 




23P2 


3980 


3996 




i'Po 




4192 




33Pi 




4300 




33P2 




4358 




l^Dl 


3820 


3800'') 


3768.9±2.5 


I'D, 


3840 


3823 




I'Ds 


3850 


3827 




2^Di 


4190 


4167") 


4159±20 


2^D2 


4210 


4195 




2'D, 


4220 


4204 





a) 



see footnote '*) to Table |T]; ^^ see footnote ") to Table |l|. 
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